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Abstract 



We study superconducting transport in homogeneous wires in the cases 
of both equilibrium and nonequilibrium quasiparticle populations, using the 
quasiclassical Green's function technique. We consider superconductors with 
arbitrary current densities and impurity concentrations ranging from the clean 
to the dirty limit. Local current conservation is guaranteed by ensuring that 
the order parameter satisfies the self-consistency equation at each point. For 
equilibrium transport, we compute the current, the order parameter ampli- 
tude, and the quasiparticle density of states as a function of the superfluid 
velocity, temperature, and disorder strength. Nonequilibrium is characterized 
by incoming quasiparticles with different chemical potentials at each end of 
the superconductor. We calculate the profiles of the electrostratic potential, 
order parameter, and effective quasiparticle gap. We find that a transport 
regime of current-induced gapless superconductivity can be achieved in clean 

superconductors, the stability of this state being enhanced by nonequilibrium. 
PACS numbers: 74.25.-q, 74.40.+k, 74.50.+r, 74.80.Fp 
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I. INTRODUCTION 



The quasiclassical theory of superconductivity was proposed by Eilenbergercl and Larkin 
and Ovchinnikov,! and later developed by Usadeli and EliashbergB It is especially useful 
to study transport in systems with characteristic length scales much greater than the Fermi 
wave length \p, which is effectively integrated out of the problem. In the seventies and eight- 
ies, the theory of nonequilibrium superconductivity was developed,! and the response of the 
quasiparticle distribution function to external perturbations was investigated. Theoretical 
developments in Refs.i@ were based on the quassiclassical Green's function (QCGF) tech- 
nique. The calculation became possible of characteristic relaxation times due to phonons 
and paramagnetic impurities when an excess of quasiparticles is injected into the super- 
conductor. In most of the cases, temperatures near the critical temperature and disorder 
strengths in the dirty limit were assumed, although the theory is potentially valid for all 
ranges of disorder and temperature. The effect of a moving condensate with a finite super- 
fluid velocity v s was studied in Refs.EHll. In particular, the quasiparticle relaxation time 
for energies between A_ and A+, where A± = |A| ± hv s pF, was calculated^ (A ± are the 
direction-dependent effective quasiparticle energy gaps; A is the conventional gap function 
playing the role of superconducting order parameter). However, the influence of v s on |A| 
through the self-consistency conditionlll was not analysed. 

During the nineties, a strong interest has developed on the physics of elastic, low current 
transport. Part of this workEl has been based on the resolution of the microscopic Bogoliubov 
- de Gennes (BdG) equations, 00 whose equivalence with QCGF was proved in the clean 
limit by Beyer et a/.0 The QCGF technique has allowed theorists to explain most of the 
experiments involving normal-superconductor (NS) interfaces. By permitting the study of 
spectral and spatial properties with inclusion of impurity averaged disorder, QCGF have 
become a most adequate tool to understand phenomena such as zero bias anomalies!! and 
re-entrance of the conduct ancejll including its non-monotonic behavior as a function of 
temperature, voltage bias, or phase difference between external superconductors.0 These 
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effects are now understood as an interplay between proximity effect and disorder.El'El 

The study of equilibrium superconducting transport for arbitrary currents was already 
initiated in the sixties.&~H Rogers^ and Bardeen0 developed a thermodynamic theory, 
and MakiS gave a unified description of the clean and dirty limits using a Green's func- 
tion method. More recently, Bagwell^! has performed a similar study for clean supercon- 
ductors within the BdG framework. The close connection between the implementation of 
self-consistency in the order parameter equation^ and current conservation has been noted 
by several authors.&il The implication is that a self-consistent description is essential in 
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transport scenarios involving large supercurrents or superconductor lengths.! 

The combined effect of a moving condensate and a nonequilibrium distribution of quasi- 
particles has been addressed recently within a self-consistent scheme.&Hil These authors 
have solved the BdG equations for one-dimensional models of ballistic transport and have 
predicted physical features such as the existence of an Andreev-transmission dominated 
transport regime,BS@ the enhancement of excess current due to a finite t> s ,00 and the 
possibility of current- induced gapless superconductivity (GS) .^SilS 

The primary purpose of this paper is to study the robustness of these effects against 
the inclusion of realistic physical factors such as scattering by impurities and the presence 
of many transverse channels. Specifically, we employ the QCGF technique to study both 
equilibrium and nonequilibrium transport in superconducting wires for arbitrary currents 
and applied voltages. Our equilibrium transport study complements work done by Makiil 
in that we compute transport properties for different disorder strengths, ranging from the 
clean to the dirty limit. In particular, we calculate the critical current density for arbitrary 
disorder, thus going beyond the Usadel equations,! valid only in the dirty limit. We also 
assume that the length of the superconductor is much smaller than the inelastic scattering 
length. 

In this work we have calculated stationary mean field solutions describing the response of 
the superconductor to an externally applied voltage. The obtained transport configurations 
are stable at sufficiently low temperatures and wire widths W not much smaller than the 



Meissner penetration length A. In this regime (specifically for W ^ A/25, according to the 



On the other hand, we assume quasi-one-dimensional superconducting wires, for which the 
Meissner effect and transverse variations of the order parameter can be neglected. This 
gives us a window of parameters for which our theory is quantitatively valid. Despite these 
considerations, we wish to point out that a satisfactory understanding is still lacking of the 
crossover from stationary configurations to dynamic phase-slips in the response to externally 
applied voltages. A unified description of both phenomena should be the object of future 
theoretical study. 

In section II, we give a brief, self-contained presentation of the QCGF technique (see 
Ref.0 for an updated review on this topic). Section III is devoted to equilibrium trans- 
port. The dependence of the order parameter, current density, and quasiparticle density 
of states on the superfluid velocity are computed exactly, and several critical magnitudes 
are calculated for different temperatures and disorder strengths. In section IV, we study 
nonequilibrium transport with elastic impurity scattering of arbitrary strength, and discuss 
the robustness of the GS regime. We solve the QCGF equation of motion with boundary 
conditions describing the injection of quasiparticles from normal reservoirs. The conclusions 
are presented in section V. 



estimate of Ref.H), both thermalU and quantumll phase-slips 



are energetically unfavored. 



II. THE QUASICLASSICAL GREEN'S FUNCTION 



We follow Refs.S0 in the derivation of the equation for QCGF. The Dyson equation for 
the matrix Green's function G in the non equilibrium Keldysh formalismEEl is 



£}G = 1 



(1) 



where 



G R G 



G 



(2) 



G A 
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Gq 1 being the usual, one-body, inverse matrix Green's function, and S the self-energy ma- 
trix in Keldysh space. Following standard notation, the symbol" indicates 2x2 matrices in 
Nambu space. The retarded, advanced, and Keldysh Green's functions are defined, respec- 
tively, as: 

G*(l, 2) = 9(tt - t 2 )[G<(l, 2) - G>(1, 2)], (3) 



G A (1,2) = -9(t 2 - ti)[G<(l,2) - G>(1,2)}, 



(4) 



G{1,2) = [G<(1,2)+G > (1,2)}, 



(5) 



where 



zG<(l,2) 



<V T (1)^}(2)> <V>t(lM(2)> 
-<^(1)V}(2)> - <^{(1)^(2)> 



(6) 



iG>(l,2) = - 



(7) 



<^ T (2)Vt(1)> <^(2)^ t (1)> 
-<V}(2)Vi(l)> -<^(2)^|(1)> 

In these expressions, we have used the standard abbreviations 1 = (ri,tx) and 2 = (r 2 ,t 2 ), 
which allow us to rewrite Eq. ([[]) as 



J d2[G \l, 2) - S(l, 2)]G(2, 3) = 5(1 - 3). 
Substracting from ([]]) its conjugate equation, we obtain 

[Gq 1 - % G\ = 0. 



(8) 



(9) 



This equation may be simplified by going to the center-of-mass and relative coordinates, 
defined as r 1)2 = R ± r/2, t 12 = T ± t/2. One Fourier transforms now with respect to the 
relative variables r and t, and introduces the QCGF defined by 



g(R,T,p,E) = - / d£G(R,T,p,E), 

7T J-oo 



(10) 



where £ = p 2 /2m—fi is the free-electron energy measured from the Fermi level and p = p/|p|. 
The key assumption in the quasiclassical approximation is that the self-energy £ is almost 
^-independent. Then we can set £ = in £ and derive for the stationary case (neglecting 
the effect of magnetic fields), 

hv F p- V R <7 = iE[f 3 ,g] - i [£,</]. (11) 

f 3 is a block- diagonal matrix with block entries like in the third Pauli matrix f 3 . After 
the subtraction procedure, a normalization condition is needed for g. In Ref.il, Shelankov 
studies the general, nonstationary case, and 

99 = 1 (12) 

is shown to be a useful choice. 

Together with the normalization condition (0), Eq. ([11]) determines the QCGF. The 



remaining physical quantities can be expressed in terms of g. Specifically, the current density 
flowing through the system is0S 

j = ~eN(Q)v F JdEj ^pTr(r 3 g), (13) 

N(0) being the single-particle density of states per spin at the Fermi level. 

Within the BCS approximation, the pairing effect is introduced in S via the order pa- 
rameter A. The resulting BCS self-energy has retarded and advanced components!! 

= -z[Re(A)f! - Im(A)f 2 ]. (14) 

In this language, the self-consistency equation for the order parameter readsll'0 



with g the electron-phonon coupling constant and Ed the usual BCS cutoff energy. Eqs. 
(|TT|), (m), (P~5|) , and ( |r5[ ) form the basic blocks for calculations in the following sections. 
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III. THE HOMOGENEOUS WIRE 



We begin by studying transport in superconducting wires with arbitrary disorder and 
homogeneous on a scale much greater than the mean free path I. It is well known that, 
up to first order in perturbation theory,!! impurities are unable to lower T c in conventional 
s-wave superconductors, the reason being that non-magnetic impurities do not break time- 
reversal symmetry.ii@ However, the effect of impurities is expected to be important when 
a non negligible superffuid velocity is present in the system. For instance, it is known that 
the critical current decreases with impurity concentration.!!! Moreover, the presence of a 
finite superflow introduces an intrinsic anisotropy in the system which one should expect 
to be sensitive to the presence of random scatterers. On the other hand, the study of 
transport in even nominally clean quasi-one- dimensional wires (I > ^ > W, £o being the 
zero temperature coherence lengthB) must include the effect of random disorder if scattering 
at the surface is diffusive.0 As mentioned above, the use of the QCGF technique permits an 
easy implementation of averages over impurity configurations. This allows for quantitative 
comparison with experiments and represents a notable advantage over techniques based on 
the resolution of the BdG equations. There have already been calculations exploring the 
extreme cases of ballistic^ and diffusive (dirty) superconductors^ in equilibrium. We wish 
to study the crossover between these two limits. In particular, we want to calculate how the 
critical current density decays with disorder and compare our results with the predictions 
based on macroscopic descriptions of the diffusive limit. 



A. The model and its solution 

We introduce disorder within the simplest approximation of incoherent multiple scat- 
tering by impurities.Ell We average over impurity configurations compatible with a given 
degree of macroscopic disorder, retaining only the leading term in an expansion in powers 
of (kpl)~ 1 ^ Under these rather standard assumptions, the contribution of disorder to the 
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self-energy may be written as 

£ imp = -—<#>, (16) 
It 

where the brackets are meant to indicate angular average over the Fermi surface [< / >= 
J" (dp/471")/]. The scattering time t coincides with the transport time r tr when the impurity 
scattering is isotropic. Eil Its precise definition is 

- = 27rn imp N(0Mp,p')\ 2 , (17) 
r 

where n- imp is the impurity concentration, and v(p,p f ) is the probability amplitude for an 
incoming electron with momentum direction p' to be scattered into direction p after collision 
with an impurity. Its square modulus is assumed to be independent of both the incoming 
and the outgoing direction. Finally, for a macroscopic description, impurity scattering is 
characterized by the disorder rate Ti/t. We will therefore use £ = Sees + Simp as the 
self-energy of our problem. 

The inclusion of a finite superfluid velocity in our set of equations leads to the addition 
of phase factors e ± * q ' R in some physical quantities. Their effect is equivalent to that of 
shifting the energy variable by an amount hvpp ■ qll@ in Eq. (O) (q is half the Cooper 
pair momentum, q = mv s /K). The presence of q leads to a nonzero supercurrent density 
and, through self-consistency equation (|15|), has a direct effect on the value of |A|. We wish 
then to solve the set of Eqs. (0), flP2D , (|T3D , and fll5D for different values of the superfluid 
velocity and K/t. 

If quasiparticles are in equilibrium with themselves and with respect to the lattice, the 
Keldysh part of the Green's function matrix ([TOj) can be expressed in terms of the retarded 
and advanced elements as g = ta,vih{E /2ksT){g R — g A ).ll These, in turn, may be written as 

g R = af 3 + (3fi 

g A = -a*t 3 + P*n, (18) 

where the scalar functions a = a(R,p,E) and (3 = /3(R,,p,E) are the generalized densities 
of states^ The normalization condition (|i2|) requires a 2 + (3 2 = 1. If one considers a bulk 



superconducting wire with periodic boundary conditions, the problem becomes uniform and 
the relevant equation of motion (ITT]) yields 



ih 

(E - hv F qu)(3 - i\A\a + —(!3<a>-a<(3>) = 0, (19) 

with u = p ■ q/q for the angular variable. This equation is solved in Appendix [A] with full 
inclusion of the w-dependence.@ Once the functions a(u,E) and (3(u,E) are obtained, we 
compute the value of |A| from Eq. (15). This new value is introduced in fll9|) to calculate 



again the generalized densities of states, from which in turn we obtain a new order parameter. 
The procedure is repeated until self-consistency in |A| is achieved. After a self-consistent 
pair potential is found for given superfluid velocity and disorder, the current density is 
calculated with Eq. ([13]). 



B. Discussion 



In Fig. [1] the order parameter amplitude |A| is plotted as a function of the superfluid 
velocity for different temperatures, and dimensionless disorder strengths T = H/tAq. One 
may note that Anderson's theoremil is satisfied at all temperatures, since |A| for v s = is 
independent of the disorder strength. For very small T, one retrieves the expected result 
that |A| ^ when v s pp/\A\ ~ 1. Increasing the disorder seems to reinforce the pair 
potential because of the much larger values of v s needed to supress A. This occurs at 
v s = Vd, with VdPF/A Q ~ 3.1, 2.6, and 1.3 for T/T c = 0.1, 0.5, and 0.9, when T = 10, 
i.e., more than twice the ballistic (small T) value. This behavior may be interpreted as the 
tendency of disorder to restore the spherical symmetry (and thus sustain the order parameter 
amplitude), counteracting the anisotropy induced by the presence of superflow (see following 
subsection) . 

This effect should not be viewed as an enhancement of superconductivity by disorder. 
Inspection of Fig. |2| shows that, as expected, the critical current density j c (defined as the 
maximum possible value of j with a nonzero A) decays with increasing disorder. This is 
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compatible with the gap behavior discussed in the previous paragraph, because the density 
of superfluid electrons is also reducedl! Fig. [| also shows that, for a given macroscopic 
current density j, there are two possible values of v s . The smaller v s yields the more stable 
configuration. Finally, we note that, for T close to T c , we reproduce the smooth behavior 
expected from Ginzburg-Landau calculations.^ 

As an illustration, we plot in Fig. ^| the values of several physical magnitudes at the point 
j = j c as a function of disorder (left pannels) and temperature (right pannels). Despite the 
above mentioned different depairing behavior, both critical j and |A| actually diminish with 
increasing disorder and temperature. However, the critical superfluid velocity is enhanced 
with T while it decreases with T. This is a manifestation of the high anisotropy needed 
to break superconductivity when disorder is strong. It also shows that disorder is not 
intrinsically depairing, while temperature is. 

Finally, in Fig. f| we plot the decay of the critical current density with disorder at low 
temperatures (T = 0.1T C ) for a wide range of T values. Inset (a) shows the linear behavior 
for small disorder (r < 0.1). Inset (b) indicates that for T > 10 the critical current follows 
a power-law (j c oc r _/3 ). Here we present a unified treatment encompassing the ballistic 
and the diffusive regimes. This requires going beyond the Usadel equations!, valid only 
within the dirty limit. Our results may be compared with the predictions of a macroscopic 
theory based on energetic arguments. If we use a phenomenological formula for the density 
of superfluid electrons (valid for I <C £o, see Ref.@), 

and equate the kinetic energy to the constant condensation energy j%/2p s , the law j c (r) = 
ic(0)/(l 4- r) 1 / 2 is obtained. This is the dotted line in Fig. |j. This simple law cannot 
reproduce the entire T dependence, and, in particular, it yields an exponent 0.5 slightly 
different from the exact (3 ~ 0.47 which we obtain numerically. 
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C. Density of states 



In Fig. H we plot the density of states (DOS) for different values of the superfluid velocity 
and the disorder strength. At low values of v s (<C v d ) we recover the characteristic BCS 
density of states. The splitting of the gap A into A + and A_ modifies the quasiparticle 
DOS. As v s raises from zero, the singularity at E = A evolves into two cusps at A±. 
The smoother character of the split singularity comes from the fact that the minimum 
quasiparticle energy depends on the momentum direction. Thus, the zero-velocity DOS 
singularity actually evolves into a distribution of singularities which, when integrated over 
the Fermi surface, yields two characteristic cusps. 

For weak and moderate disorder (r = 0.1 and 1), the GS regime can be achieved with 
values of v s within the stable branch (see Fig. for the corresponding temperature T/T c = 
0.1). This result is important, since it shows that the GS state occurs in a stable manner 
for equilibrium transport in relatively clean superconducting wires. GS can also be reached 
in dirty (r = 10) superconductors but only for unstable values of v s . 

Comparison of the DOS curves for T = 1 and T = 10 in Fig. [| shows that, for some 
values of v s pf/A (such as 1 or 1.5) for which one would expect to have GS, the effect 
of disorder is that of restoring the gap. This effect may be understood as resulting from 
the directional randomization induced by multiple impurity scattering. The idea that GS 
should exist for v s > \A\/pp comes from a kinematic analysis that applies to plane waves 
or to quasiparticle states differing little from them. This is the case of weak disorder. By 
contrast, in strongly disordered superconductors, the exact quasiparticle states necessarily 
involve a strong mixture of plane waves pointing in many directions. Disorder helps to 
preserve the gap in the exact DOS because, at very low energies, semiclassical quasiparticle 
trajectories can only select plane waves from a narrow solid angle of gapless directions. Not 
being able to mix many momentum directions, quasiparticle states at such low energies cease 
to exist. This effect is analogous to the appearance of a minigap in the DOS of a diffusive 
normal metal in contact with a superconductor.0'@ 



11 



IV. NONEQUILIBRIUM TRANSPORT IN A DISORDERED WIRE 



The last section of this paper is devoted to nonequilibrium transport in a disordered 
superconductor symmetrically connected to two normal reservoirs with different chemical 
potentials, as schematically depicted in Fig. |6|. This boundary condition introduces a 
nontrivial spatial dependence of the physical magnitudes and makes the equation of motion 
([TT|) more difficult to solve. In particular, the Keldysh component of the Green's function 
matrix (IT]) contains now all the relevant information on the quasiparticle local distribution 
function and on the momentum-relaxation processes taking place within the superconductor 
due to disorder. Of course, the main new physical ingredient involved in our calculation is 
the presence of a phase gradient in the order parameter, which has to be determined self- 
consistently at each point for every value of the chemical potentials at the reservoirs. Like 
in Refill, the loss of translational invariance makes it harder to compute the pair potential, 
which now has to be determined both in its real and imaginary parts, as demanded by a 
locally self-consistent calculation. 

We are interested in the behavior of the new transport regimes which were mentioned 
in the Introduction. In Refs.BS it was shown how the regime dominated by AT requires 
the presence of moderately reflecting barriers at the NS contacts in order to be realized. 
Unfortunately, the mathematical implementation of the appropriate boundary conditions 
in the context of the QCGF is very cumbersome, since it involves the solution of non- 
linear equations.0 Due to this complication, we have chosen to perform the study for ideal 
NS contacts (where the boundary conditions reduce to continuity of the QCGF at the 
interfaces) , leaving for the future a more systematic calculation for arbitrary barriers at the 
interfaces, in the spirit of Ref.0. With this assumption, the AT regime cannot be studied, 
and we will focus instead on the sensitivity of the GS regime to the degree of disorder in the 
superconducting wire. 
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A. Solution of the kinetic problem 



We assume a superconducting wire of length L, connected at its ends z = ±L/2 with 
perfect normal leads (characterized by having A = 0). These normal leads are connected to 
large normal reservoirs through ideal contacts, so that the chemical potentials characterizing 
the population of incoming electrons and holes are those of the reservoirs from which they 
were injected.il Periodic boundary conditions are considered in the transverse directions, 
perpendicular to the transport direction z. For a given voltage V = (fi^ — /ijj)/e, we wish 
to solve again Eq. ([11]) with the normalization prescription ([12]) for a specific shape of the 
order parameter A(z) = \A(z)\e lip ( z \ Once this is done, the solution is used to compute 
a new A(z), and the procedure is repeated until self-consistency is achieved. Finally, the 
current density j is calculated from Eq. (|l"3|) . 

In the stationary limit, and for an arbitrary shape of A, it is possible to parametrize the 
retarded and advanced QCGF asll 

g R = af 3 + P?! + 7f 2 

g A = -a*f 3 + /rf 1 + 7 *f 2 , (21) 

where a, j3, and 7, are now scalar functions of the longitudinal coordinate z (as well as 
of u and E). When Eq. ( pT[) is introduced into the equation of motion ( |TT]) a new set of 
equations is obtained: 

dct fi Ti 
huv F — = (i- <7> -2iImA)/3 - (i- </3 > +2iReA)7 

OZ T T 

88 h h 

huv F ^- = (-i- < 7 > +2ilmA)a + (i-<a> +2E)^y 
oz t r 

O'y ti Ti 
huvp- 1 - = (i-<p>+2iReA)a - (i-<a>+2E)p, (22) 

OZ T T 



where, as usual, the brackets stand for angular average. The normalization condition ( |T2"D 
implies a 2 + 3 2 + 7 2 = 1. In Appendix [B] we solve this set of equations with the appropriate 
boundary conditions. 

While the retarded and advanced parts of the equation of motion still give us the general- 
ized densities of states of the problem, the Keldysh part g of (|l0f), which contains information 
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on actual occupations, is no longer trivial. Being in a nonequilibrium context, we need to 
write its corresponding equation of motion, which reads 

Huv/l = iE(f 3 y - M - + ttf - i«± - (23) 

One way of parametrising g which automatically satisfies Eq. (|12|) is 

g = g R h - hg A , (24) 

with h an arbitrary distribution matrix. Schmid and Schoni proposed a distribution matrix h 
of the diagonal form h = (1 — 2fi) — 2/rf 3 . This procedure is extremely useful in situations 
where the dirty limit is valid and the boundary conditions are independent of u, but 
it does not permit a complete separation of the densities of states and the quasiparticle 
distribution function.il An alternative parametrization was suggested by Shelankov,i 

9 = \ {fiA R Pi A + / 2 i^A) , (25) 

where /x,2 are scalar functions of the variables (z, u, E), and P^ 2 = 1 ^9 R an d P\i = 1 Z F9 A 
satisfy 

pR(A)pR(A) _ 2 pR(A) . _ 1 2 

pR(A)pR(A) = Q (26) 

Multiplying Eq. ( p^) by P^P^ on the left, taking the trace, and using the cyclic properties 
as well as Eq. (pE|), one may find after some algebra, 

huvFTriPfP*)^- = -iTr{PfP*[(Z R - t A )j\ - S]}. (27) 

Eq. ( pTj) is a closed expression that permits to obtain f±(z, u, E) since, due to the symmetry 
of the problem, f 2 (z,u,E) = —f\{—z,u,E). Actually, (1 — /i)/2 coincides with the usual 
Fermi distribution function in equilibrium, and with the quasiparticle distribution function 
within a semiconductor model.il To fix the population of incoming electrons from each 
reservoir, it is natural to use the boundary conditions 
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fi(-L/2,u,E) = tanh[(£ - eV/2)/2k B T\, ifw>0 
fi(L/2,u,E) =ttmh[(E + eV/2)/2k B T], ifu<0. (28) 

The resolution of Eq. (B^) with this boundary conditions also requires a self- consistent 
determination o/S. The contribution of Ei mp has to be calculated at each energy from the 
solutions of Eqs. ([22]) and ( p7|) . The corresponding angular average of the whole g is then 
performed and introduced in Eq. ([16]) until one achieves a self-consistent value of £i mp . 
Once this is done at every energy E, the self-consistent procedure for A may be initiated. 

Finally, when the self-consistent solutions for a given voltage are found, one may calculate 
the electrostatic potential in the structure asl§ 

e<j>{z) = ~JdE<Trg>, (29) 

which is directly related to the electronic density. Eqs. ( f29|) and fll3|) will be used to 
compute the electrostatic potential profile and the current density which we discuss in the 
next subsection. 

B. Discussion 

In Fig. [7] we plot the variation of the current density with the applied bias V between the 
normal reservoirs for different values of the disorder. The variation of the order parameter 
amplitude at the center of the NSN structure z = is represented in the lower figure. 
Throughout this subsection, we use the values of L = 3£ for the length of the superconductor 
and T = 0.01T C for the temperature. One may see how when disorder is present (r ^ 0), 
there exist two regimes with different current-voltage slopes. For low bias, dj/dV is always 
bigger than its corresponding normal value (which is that attained at higher voltages). This 
fact shows clearly the different effect of diffusive (dirty) regions and tunnel junctions on 
superconducting transport. In a purely ballistic conductor attached via tunnel barriers to 
the reservoirs, the conductance would always be bigger in the normal statejl the main reason 
being the much smaller probability for the simultaneous tunneling of two electrons forming 
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a Cooper pair. However, if the contacts are good enough, and scattering is dominated by 
spatially distributed impurities, the superconducting state supports a greater amount of 
current due to the lack of normal reflection of quasiparticles decaying into Cooper pairs. In 
a perfectly clean sample, the superconducting and the normal states cannot be distinguished 
in what regards to transport (both display the same slope in the T = case).ili3@@ 

There is a correlation between the slope discontinuity for the current and the vanishing 
of A (z = 0). This is particularly clear for r > 1 since then the effect is more marked 
and occurs at lower voltages (eV ~ 1.3Ao). This is the point at which the whole structure 
becomes normal and transport is governed by the Boltzmann equation. The case of weak 
disorder (r = 0.2) is a bit different. One may notice that the current density goes down 
slightly above eV = 1.5Ao, while A(0) remains finite for still higher voltages. The structure 
is still superconducting (in the sense of having a nonzero order parameter) but the current 
is practically indistinguishable from that of the normal state. The superconductor is in the 
GS regzmi^SiH (see discussion below). 

In Fig. |3| one may compare in more detail the differences between the situations with 
weak and strong disorder. We will concentrate on the different behavior of the GS regime. 
Several physical magnitudes have been represented as a function of position for different 
values of the applied voltage. The top pannels represent the electrostatic potential along 
the superconductor S. Due to the interfaces with normal reservoirs and to the presence of an 
applied voltage, there is a penetration of the electric field within S, causing the quasiparticles 
to have a chemical potential different from that of the condensate (equal to zero in these 
graphs) .0 

The relaxation of the electric field within a superconductor is one of the most heavily 
studied topics in the literature on nonequilibrium superconductivity.il The physics we en- 
counter here is somewhat different because we deal with a finite superconductor without 
inelastic scattering, connected symmetrically to normal reservoirs. Two different elastic 
mechanisms contribute to the spatial relaxation, after length A, of the excess charge density 
due to quasiparticles near the boundaries. One may have normal reflection of quasiparticles 
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returning to the reservoir. This translates into a law [1 + {z + L/2)/l}~ 1 for <p(z) on the 
left boundary, characteristic of decay by incoherent multiple normal scattering.ElH On the 
other hand, beyond a certain length, quasiparticles cannot propagate because their energy 
is insufficient to overcome the effective gap, since this increases as one gets deeper into the 
superconductor. This occurs at a characteristic length from the boundary hvp/eiV), where 
e(V) is an effective energy that diminishes with V because of the self-consistent decrease 
of |A| in all the structure. This is the length scale at which quasiparticle conversion into 
Cooper pairs by Andreev reflection takes place. Andreev scattering generates a quasiparticle 
of opposite charge and thus tends to quickly suppress the charge excess.0 These considera- 
tions give a good account of the results for <p(z) in Fig. which shows a roughly constant A 
as a function of the applied voltage up to the transition to the normal state for T = 5 (dirty 
limit, I dominates), and a continuously increasing A when T = 0.2 (clean limit, hvp/eiV) 
dominates). When the structure becomes normal, the voltage profile is that which results 
from substracting the curves [1 + (L/2 ± z)/V\~ l . In the very clean limit (I 3> L) this results 
in an essentially straight line. 

Center and bottom pannels represent the |A| and A_ profiles, respectively. Since tp'(z) 
is no longer constant, A_(z) = |A(z)| — Hvptp'(z)/2. When A_ < at some point, the 
superconductor becomes locally gapless, at least for plane waves.il In particular, near the 
interfaces, A_ becomes negative very quickly as V increases, due to the smallness of |A|, 
while it remains positive in the bulk (near z = 0). As the voltage increases, one may 
note a different evolution of the left and right graphs before the normal state is reached 
(this occurs at eV = 1.86A and eV = 1.35A , respectively). In the diffusive regime, the 
superconductor cannot be globally gapless. All the A_ curves are positive in the central 
region. Quasiparticles with energies below the effective gap penetrate the superconductor 
partially and become Cooper pairs. However, if A_(0) < 0, this transfer process becomes 
marginal. Normal scattering dominates and, due to the strong disorder, the quasiparticles 
can hardly be transmitted across the system. The directional randomization^ results in a 
sharp transition to the normal state. On the contrary, for T = 0.2, Fig. || shows that it 
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is still possible to have A_ < for all z (for, e.g., eV = 1.5Aq), while the system is still 
superconducting. The quasiparticle density of states is nonzero at all energies because the 
weak disorder is unable to randomize the quasiparticle nonequilibrium distribution, which 



on the wire length and on its effective dimensionality^) , the GS becomes unstable and the 
system goes normal. 

Since realistically clean superconducting wires can be obtained nowadays, it could be 
possible to measure finite superflow effects passing a large current along a given direction of 
a superconducting sample, while probing the gap and/or the density of states with a weak 
tunnelling current. 



We have used the technique of the quasiclassical Green's function (QCGF) to study the 
effect of disorder and many transport channels in situations where superconductors present a 
non-negligible superfluid velocity. The QCGF technique has allowed us to include impurity 
scattering in a tractable manner, and to perform realistic calculations of disorder-averaged 
physical quantities. 

The first part of this article has been devoted to the study of equilibrium transport in 
homogeneous wires. We have seen that disorder tends to restore the spherical symmetry 
and strengthen the order parameter amplitude, in contrast with the pair breaking effect of 
finite condensate flow. However, its net effect is that of diminishing the current. The critical 
current density has been calculated for disorder strengths ranging from the clean to the dirty 
limit. The self-consistent density of states for different values of the superfluid velocity and 
disorder has also been discussed. We find that, due to its direction mixing effect, sufficient 
disorder restores the quasiparticle gap in transport contexts where one would expect gapless 
superconductivity (GS). 

The second part has focussed on nonequilibrium transport in normal-superconductor- 



results in the preservation of superconductivity^ For sufficiently high voltages (depending 





V. CONCLUSIONS 
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normal structures. Clean contacts between the superconductor and the normal reservoirs 
have been assumed, with boundary conditions that describe incoming quasiparticles from 
reservoirs with different chemical potentials. We have found that the current-induced GS 
regime is very sensitive to the presence of disorder. It is suppressed in the diffusive limit, but 
is stable in sufficiently clean superconductors. We have calculated the spatial profiles of the 
electrostatic potential, order parameter amplitude, and energy threshold for quasiparticle 
transmission. We have found that they show subtle differences in the ballistic and diffusive 
regimes in what regards to the current- induced phase transition to the normal state. 
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APPENDIX A: IMPURITY- AVERAGED GENERALIZED DENSITIES OF 

STATES 

Eq. (^) is actually an algebraic equation if one considers the averaged c(E) =< a(u, E > 
and s(E) =</3(u, E) >as known quantities. If it is so, one may define t(u, E) = /3/a, whose 
formal solution is 

V ' ' (E -hv F qu + ihc{E)/2r) K ' 

Now, c(E) and s(E) may be found if one notices the normalization prescription a 2 + /3 2 = 
1, which makes a = (1 + t 2 ) _1//2 and (3 = t(l +t 2 )~ 1 / 2 . If we change variables in the angular 
integrals for c(E) and s(E): 

<•*- i(|A| +hs(E)/2T) dt 



c(E) 



2hv F q t 2 (l + t 2 )V2 
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s(E) 



i(\A\ + hs(E)/2r) dt 



(A2) 



/*+ 2hv F q t(l + t 2 ) 1 /2' 

with t±(E) = i(\A\ + hs(E)/2r)/(E ± hv F q + ihc(E)/2r). Integrals in Eq. fl£2|) give two 
coupled equations for c(E) and s(E) 



i(\A\+hs/2r) 
2TiVpq 



1+^)1/2 (1 + ^)1/2 



i(|A|+ns/2r). t-fl + a + d) 1 / 2 ] 

S = — — ! 1 log 



(A3) 



2hv F q "°t+[l + (1 +tl)V2]' 
where we have omitted the dependence on E for simplicity. These equations can be solved 
self-consistently at each energy using standard techniques. 



APPENDIX B: THE SCHOPOHL-MAKI TRANSFORMATION 



Eqs. 



can 



be decoupled using the Schopohl-Maki transformationliE 



2/1,2 



1 + a 



(Bl) 



which leads, with the help of Eq. (12), to the following Riccati differential equations: 



huv F ^- - 2iEy x + A 2 yf - A x = 
oz 



huv F ^ + 2%Ey 2 - A x y\ + A 2 

oz 



0. 



(B2) 



where E = E + ih<a> /2r, A x = A + h <(3 - 27 > /2r, and A 2 = A* + h</3 + ry> /2r. 
When one solves Eq. (p2|) in the normal leads (where A = and T = 0), the solutions 



arc 



yi = yi (z )e 2 ^ z - z ^ huVF 

y2 = y2(z )e~ 2 * Eiz - z °y^, 



(B3) 



with z Q = ±L/2. In the imaginary time representation, one has the freedom of choosing 
either the positive or the negative imaginary axis for the Matsubara energies involved in the 
problem. Considering the usual criterion of E = iu n , with uj n > 0, one obtains as appropriate 
boundary conditions that ensure the finiteness of the solutions in the normal leadsS 
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Vi{-L/2) = y 2 (L/2) = if«>0 

y 2 {-L/2) = yi {L/2) = ifu<0. (B4) 

On the other hand, the symmetry relation A(z) = A*(— z) guarantees the following sym- 
metries for the solutions of the generalized densities of states (regardless of the voltage 
difference between the normal reservoirs): 

a(z, u, E) = a(—z, u, E) 
(3{z,u,E) = f3{-z,u,E) 

1 (z,u,E) = - 1 (-z,u,E), (B5) 

which translates into 

Vi(z,u,E) =y 2 (-z,u,E). (B6) 

Eq. (|B2|) reduces then to one effective equation for each y, which can be solved numerically, 
and the generalized densities of states may finally be calculated from the relations 

l-Vi(z)Vi(-z) 



a(z) 
(3(z) 



i + yi(z)yi{-z) 
yi(z) +yi(-z) 

1 + Vi(z)vi(-z) 



1{Z) = y - "■(--)] , (B7) 
l+yi(z)yi(-z) 

where we have omitted the u and E dependence. As explained at the end of subsection 
IV.A, a self-consistent procedure is also needed at each energy to obtain the right angular 
averages entering the equations of motion. 
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FIG. 1. Impurity averaged order parameter amplitude as a function of the superfluid velocity, 
for different values of the normalized disorder rate T = H/t Aq and temperature. 
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FIG. 2. Same as Fig. |l] for the current density. 
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FIG. 3. Impurity averaged critical magnitudes as a function of the disorder (left) and of the 
the temperature (right). 
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FIG. 4. Impurity averaged critical current density as a function of the disorder for T = 0.1T C . 
Insets: (a) Ballistic case (r < 1): Linear decay of the critical current with the disorder rate, (b) 
Diffusive case (T » 1): Power-law j c oc r _/3 . We find /3 ss 0.47, slightly different from the simple 
prediction of 0.5 (dotted line). 
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FIG. 5. Impurity averaged density of states for several superfiuid velocities and disorder 
strengths. All results are given for T = 0.1T C . 
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FIG. 6. Schematic representation of a typical NSN structure, where S is a dirty superconduc- 
tor. The normal leads are assumed to be perfect and connected to wide reservoirs through ideal 
contacts. Grey colour marks the superconducting zone of length L. Small crosses signal possible 
impurity positions. The picture shows a specific realization of disorder within the superconductor. 
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FIG. 7. Top: Current density as a function of the applied voltage in the structure of Fig. 
||, for different values of the disorder. Bottom: Order parameter amplitude at the center of the 
superconductor, for the same values of T. Here, L = 3£o and T = 0.01T C . 
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FIG. 8. Comparison between ballistic (left) and diffusive (right) limits in the structure of Fig. 
^. Top pannels: Profile of the electrostatic potential within the superconductor. Center pannels: 
Order parameter amplitude. Bottom pannels: Threshold A_ = |A| — Tivpip'j'l = |A| — mvFV s . At 
positions where A_ < the superconductor is gapless.H L and T as in Fig. ^ 
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